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First principles density functional theory (DFT) is used to investigate the electronic structure 
of /3-Mn02- From collinear spin polarized calculations we find that DFT+Ueh predicts a gapless 
ferromagnet in contrast with experiment which indicates an insulating antiferromagnet. The inclu- 
sion of anisotropic Coulomb and exchange interactions in the DFT+U approach, defining U and J 
explicitly, corrects these errors and leads to an antiferromagnetic ground state with a fundamental 
gap of 0.8 eV consistent with low temperature experiments. To our knowledge, this work on /3-Mn02 
represents the first demonstration of a case in which the application of fully anisotropic interactions 
in DFT+U determines the magnetic order and consequent band gap, while the more commonly used 
effective U approach fails. Such effects are argued to be of importance in many insulating materials. 
The mechanism leading to an increase in band gap due to anisotropic interactions is highlighted 
by analytical calculation of DFT+U d-orbital eigenvalues obtained within a Kanamori-type model. 
Magnetic coupling constants obtained by the fitting of a Heisenberg spin Hamiltonian to the energies 
of a range of magnetic states assist in rationalizing the finding that anisotropic interactions enhance 
the stability of the experimentally observed helical antiferromagnetic order. The plane wave PAW 
method yields poorer results for the exchange couplings than full-potential LAPW calculations. Fi- 
nally, we compare the DFT+U results with exchange couplings obtained from hybrid functionals. 
It is argued that anisotropic interactions should be included in DFT+U if the results are to be 
properly compared with those from hybrid functionals. 

PACS numbers: 71.27+a, 71.15.Mb, 75.50.Ee 



I. INTRODUCTION 



Density functional theory has been used with success 
in understanding the properties of many materials. How- 
ever, in the class of transition metal oxides the perfor- 
mance of traditional DFT methods, typically based upon 
local spin density or generalized gradient (GGA) approx- 
imations, is often found to be inadequate. For instance, 
these methods predict metallic states for many insulating 
transition-metal oxides^. This is understood to be due 
to the presence of strong on-site Hubbard U interactions 
resulting from the electronic correlations in open d and 
/ shells. The DFT+U approach has been introduced^ 
specifically to overcome such deficiencies. This approach 
incorporates an enhanced treatment of Coulomb and ex- 
change interactions in self-consistent DFT, improving the 
description of such materials in many cases. The present 
study of /3-Mn02 makes the further point that an accu- 
rate treatment of fully anisotropic Coulomb and exchange 
interactions may be of vital importance in application of 
the DFT+U method, proceeding beyond the conceptu- 
ally simpler 'effective U' approach that is often applied. 

The material /3-Mn02 is of technological and fun- 
damental scientific interest, being the most stable and 
abundant polymorph of the M11O2 oxide family. Techno- 
logically it is important as a cathode material in lithium 
ion batteries^ and as a catalyst^. Experimentally (3- 
Mn02 occurs in the rutile structure^ with space group 
PA/mnm (#136) and lattice parameters a = 3.3983 
A and b = 2.8730 A. It orders magnetically at TV = 



92 K-i. Below this temperature, transport measurements 
indicate an enhancement of the resistivity up to large 
values (~ 10 5 fi cm) at T =0 K, indicative of insulat- 
ing behavior—. The magnetic transition has been investi- 
gated in depthSr— and results in a helical antiferromag- 
netic order with pitch (7/2)c. For simplicity, this pa- 
per initially focusscs on collinear antiferromagnetic order 
between neighboring Mn sites along (111). This corre- 
sponds to opposing spins on the center and corner man- 
ganese sites in the unit cell shown in Fig.QJa). Antiferro- 
magnetic coupling of these sites has been shown to be key 
to the helical magnetic orderi^. Following this we return 
to treat the non-collinearity explicitly. In Fig. QJb) the 
helical order of pitch (7/2)c in the magnetic unit cell is 
shown. 

Two DFT+U approaches are compared in this work. 
The first we term the 'full anisotropy' approach, where 
U and J are both defined distinctly and the resulting 
Coulomb and exchange matrices respect the full orbital 
dependence of the d-d interactions. The second, intro- 
duced by Dudarev et aiJ^, is denoted the 'effective U' ap- 
proach where only U is distinctly defined as Ues = U — J , 
and J is set to zero throughout. 

It is important to note that a large majority of DFT+U 
calculations in the literature are performed with the effec- 
tive U and therefore this work is of importance to many 
researchers. In some codes only the effective U approach 
is implemented (CASTEP— ), while in others it is the 
default (VASP— ) or is recommended (Wien2ki£). 

Previous electronic structure calculations of /3-Mn02 
have been performed by Franchini et al£. Similar to the 
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FIG. 1. (Color online) (a) The structural unit cell of /3-Mn0 2 
showing the approximate octahedron formed by the oxygen 
(small red spheres) coordination of the central manganese 
(large purple spheres). First four nearest neighbor magnetic 
couplings also indicated, (b) The magnetic unit cell with 7/2 
structural units along the c-axis. The magnetization vector 
is illustrated for the helical order from Yoshimori^. Note the 
atoms are included to indicate scale, but do not represent an 
integer number of structural unit cells. 



tonian employed was adapted to a cubic symmetry that 
/3-Mn02 does not possess. Our work is the first to study 
the energetic stability of the full non-collincar order. 

Previous studies have demonstrated that anisotropic 
interactions in DFT+U may give rise to spin ordering of 
orbitals^ and influence the canting angle in non-collinear 
magnetic systems^. 



II. THE DFT+U CORRECTION 

The DFT+U method in essence introduces a correction 
to traditional DFT energies based upon the occupation of 
correlated orbitals. Typically these orbitals correspond 
to strongly localized d and / states, wherein the on-site 
Coulombic and exchange interactions between electrons, 
usually described by U and J, are not adequately treated 
by traditional functionals. We may write this correction 
to the total energy in the form: 

AE = Ej - E dc (1) 

The second term is a double counting correction which 
has been discussed in detail previously^, while the first 
corresponds to the interactions between the correlated 
electrons: 



present work, they found that the effective U method pre- 
dicts a metallic ferromagnetic (FM) ground state over a 
wide range of U values. This contravenes the experi- 
mental evidence^ which is, in turn, in better accord with 
the insulating antifcrromagnctic (AFM) ground state ob- 
tained by Franchini et al. using hybrid functionals. It is 
noted, though, that the latter study did not consider the 
non-collincar magnetic order, but did suggest that the 
deficiency of the effective U approach may be due to the 
fact that the improved treatment of correlations is lim- 
ited to localized d-states. In this study we show quan- 
titatively that the inclusion of anisotropic Coulomb and 
exchange interactions in DFT+U corrects the effective U 
description. 

An early periodic Hartree-Fock (HF) study by Mack- 
rodt and Williamson^ found a very large gap of ap- 
proximately 4 eV, the FM ground state being favored 
over AFM by 60 meV per formula unit. Other stud- 
ies by Balachrandan et al. have also considered /3- 
Mn02 as part of wider studies relating to polymorphic 
Mn02 as a lithium battery insertion cathod o 17 i 18 . A re- 
cent DFT+dynamical mean field theory (DFT+DMFT) 
studyi^ also found a gap in the excitation spectra of AFM 
/3-Mn02 within the fixed experimental structure. How- 
ever, the total energies of likely magnetic states were not 
compared. Furthermore the magnetic order was treated 
within an "appropriate cubic" cell that does not fully rep- 
resent the experimental magnetic order, while the Hamil- 



Ei = \ J2 W ™n>n m *n ml<T , (2) 

ma+m'cr' 

Here, for simplicity and clarity we write the sum in 
the Hilbcrt space in which the occupations are diago- 
nal, which can always be achieved by a rotation. The 
Coulomb matrix W™ m , may be expressed as: 

^ T mm' — U mrn ' J mm' ^ oo' (3) 

Umm' are the direct Coulomb interactions which oper- 
ate between states irrespective of spin and J mm' & a a' de- 
fines the spin-dependent exchange interactions. The ele- 
ments of the two interaction matrices are all derived from 
the parameters U and J taken as input to the DFT+U 
method; however, what is often not appreciated is that 
the direct interaction U mm > is determined not only by U, 
but also by J. 

It is J that determines the anisotropy of the Coulomb 
interactions^ in the DFT+U scheme, via a linear com- 
bination of Slater integrals. In essence these interactions 
have a tendency to increase the mutual distance of each 
pair of electrons and therefore promote Hund's second 
rule. Approaches to DFT+U based on an effective U 
omit this anisotropy by setting J to zero resulting in 
Umm' = U and J m m' = for m ^ m! . In contrast, with 
full anisotropy the Coulomb and exchange matrices (for 
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U = 6.7 and J = 1.2 eV) arc: 
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The off-diagonal elements of the Coulomb U m m' and 
exchange Jmm 1 matrices vary significantly due to the 
anisotropy of the interactions within the localized d- 
manifold. 

The input parameters U and J are translated into ma- 
trices Uram 1 and Jmm' by appealing to the atomic physics 
of the ion under HF theory^. This well-established 
method provides an accurate description of interactions 
in isolated atoms2£ (excepting those processes that in- 
volve a change in the number or significant change in the 
character of pair correlations, e.g. spin excitations) and 
is therefore modelled within the DFT+U approach. HF 
methods, including on-site hybrids^, inherently involve 
the computation of integrals over the occupied wavefunc- 
tions to determine the Slater integrals F°, F 2 , F 4 and 
F 6 . A linear combination of such integrals then deter- 
mines the Coulomb and exchange matrices. However, in 
the DFT+U method, rather than calculating the Slater 
integrals explicitly we parameterize them using U and J. 
Thus, for the d-orbital manifold of interest here, F Q = U 
and, guided by atomic calculations, we set F 2 +F 4 = 14 J 
with the condition that F 2 /F i = 8/5. Higher order inte- 
grals occur only for /-electrons. One advantage of such 
a parameterization is that we may directly apply values 
of U and J that are more representative of the effects of 
screening due to both correlations and interaction with 
the crystalline environment, with renormalized U values 
typically in the range 4-8 eV, whereas a direct HF cal- 
culation of U = F° provides an unrealistic value of ap- 
proximately 22-26 eV. 



III. ELECTRONIC STRUCTURE RESULTS 

We have calculated the electronic structure using 
GGA+U. The VASRii code was employed using the 
PAW method with a plane- wave cut-off of 600 eV. A min- 
imum 6x6x6 k-point grid was used in the full Brillouin 
zone. For comparison, key results were also computed 
separately using the full potential LAPW approach in 
Wien2ki£ with RK ma x = 7.0 and muffin tin radii of 2.01 
ao for Mn and 1.51 ao for O. The fully localized limit 
(FLL)22 double counting correction was employed which 
is considered approprite for localized d-states in /3-Mn02- 

We directly obtain Ues by the ab initio method as ap- 
plied in Wien2k2£, yielding a screened value Ues = 5.5 
eV as might appropriately be applied within the spher- 
ically averaged effective U approach. Since the Slater 
integrals, F 2 and F 4 , are typically only weakly screened 



TABLE I. Unit cell dimensions, unpaired spin moments and 
energies of antiferromagnetic (AFM) and ferromagnetic (FM) 
ordering for /3-Mn02 calculated with VASP. Experimental 
structural data from Regulski el al£&. Total energy for ferro- 
magnetic order is set as the zero reference. 
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4.496 


4.443 


4.451 


4.398 


C (A) 


2.955 


2.991 


2.936 


2.940 


2.874 


H (p B ) 


3.24 


3.61 
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in solidsi"^, we use the atomic-limit^ value for Mn 4+ 
J = 1.2 eV. Therefore, in the full anisotropy GGA+U 
approach, we set U = 6.7 eV and J = 1.2 eV. In work to 
be published elsewhere, we show that these values yield 
accurate lithium intercalation voltages for /3-Mn02 act- 
ing as a cathode in lithium ion cells. 

In Table U we compare the relative total energies ob- 
tained for FM versus AFM ordering from VASP. The 
effective U approach gives the lowest total energy for a 
FM spin configuration with a stabilization of 48 meV 
per formula unit. In contrast, full anisotropy GGA+U 
shows an energetic stabilization of 20 meV per formula 
unit for AFM order. A qualitatively similar result is also 
obtained using the local spin density approximation + U 
(LSDA+U) method, indicating that the relevant physics 
is not specific to a particular choice of local or semi-local 
functional. We have verified in Wien2k using GGA+U 
that AFM is stable (by 35 meV per formula unit) only 
when full anisotropy is included, leading to the conclu- 
sion that this effect is not implementation or basis set 
dependent. However, applying the effective U approach 
with the LAPW method in Wicn2k leads to a magnetic 
energy difference of just 2 meV per formula unit; the two 
magnetic orders are essentially degenerate. We return to 
discuss this in section UlI Bl where the exchange constants 
dictating the magnetic ground state are calculated. 

We have also tested the AFM stabilization energy for 
other reasonable ranges of Ues (4-6 eV) and J (0.9 
- 1.4 eV) and find that, while the magnitude varies, 
AFM order is always favored by application of the fully 
anisotropic approach. Most importantly, the large differ- 
ence in the relative energies of FM and AFM states is 
robust. 

As is common with the GGA+U method, Table U 
shows that the cell parameters are overestimated, but the 
agreement with experiment is improved when anisotropy 
and the lower energy AFM order are both employed. 
The unpaired electron count of 2.95 obtained with full 
anisotropy is consistent with the unpaired spin density of 
a Mn 4+ = d 3 =>■ 3/j.b valence state formally expected for 
M11O2. As discussed by Ylvisaker— et at the inclusion 
of J can lead to such lower energies for smaller magnetic 
moments by promoting spin and orbital polarization. We 
note that the calculated unpaired electron count of 2.95 
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FIG. 2. (Color online) Density of states of /3-Mn02 from 
Wien2k. Panel (a) effective U ferromagnetic, (b) is full 
anisotropy ferromagnetic and (c) full anisotropy antiferromag- 
netic. (f.u.=formula unit.). 



does overestimate the experimental^ low temperature 
ordered moment of 2A^b- These experiments include 
contributions from orbital moments and spin excitations 
that are not captured by our calculations. Furthermore, 
there is a large phase space for spin fluctuations avail- 
able in this system due to the large number of competing 
magnetic orders*^. To capture the effects of these fluc- 
tuations a theory with dynamical frequency dependent 
interactions may be required^. 

In Fig. [2] we show the electronic density-of-states 
(DOS) computed for /S-MnC^. The top two panels pro- 
vide the DOS in FM order, Fig. dja) showing the re- 
sult for the effective U approach while Fig. [2fb) includes 
anisostropy in the GGA+U treatment. The anisotropic 
J interactions introduce a small gap while the effective 
U approach yields a metallic density-of-states. Fig. [3Jc) 
shows the DOS for the AFM spin configuration obtained 
from GGA+U with full anisotropy, wherein a band gap 
of approximately 0.8 eV in width is evident, consistent 
with experimental measurements^. 

It is instructive to see what the effect on the underlying 
spin and charge density is due to anisotropic interactions. 
The Mn local environment has two apical Mn-0 bonds 
that are at 90° to the four in-plane Mn-0 bonds together 
comprising the approximate octahedron. However, the 
angle between the in-plane oxygens deviates from 90° 
to values of approximately 100° and 80° for consecutive 
pairs as shown in the inset to Fig. [3] Consequently, the 
octahedral symmetry is approximate and to obtain the 
orbital projected partial charges we set the z-axis along 
the apical Mn-0 bonds and the x,y-axes within ~ 5° 
from the in-plane Mn-0 bonds. In Table [n] the partial 
charge and net spin for the Mn d-orbitals is shown. The 
results indicate that the crystal field dominates the or- 
bital occupations giving a clear e g to t2 g splitting and 
Hund's first rule is obeyed with majority spin t 2g occu- 
pations close to unity. 



TABLE II. Projected orbital charge and net spin occupations 
on the Mn site for an antiferromagnetic spin configuration of 
/3-Mn0 2 obtained with Wien2k. 

Effective U Full Anisotropy 
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Spin 


Charge 


Spin 


d z 2 


0.609 


0.199 


0.599 


0.090 


d x 2 — y 2 


0.640 


0.183 


0.634 


0.075 


dxy 


0.965 


0.875 


0.965 


0.879 


dxz 


1.021 


0.839 


1.020 


0.843 


dyz 


1.008 


0.849 


1.007 


0.854 



The partial charges in Table [TT] obtained with full 
anisotropy and effective U are very similar, all deviations 
being less than 0.01. In contrast the net spin per orbital 
varies significantly between the two approaches. For ex- 
ample the d z 2 net spin is 0.199 for effective U and 0.090 
for full anisotropy. The deviations are also large for the 
other e g orbital, d x 2_ y 2 , with a net spin of 0.183 for effec- 
tive U and 0.075 for full anisotropy. However, for the t^g 
orbitals d xy , d xz and d yz there is little difference in the 
spin occupations between effective U and full anisotropy. 
It is the net spin of the e g orbitals that is reduced by 
full anisotropy. Small changes to the spin polarization 
of these predominantly bonding orbitals can have a large 
impact on the magnetic exchange. The full anisotropy 
approach, by respecting the full symmetry of the d-d in- 
teractions gives a better description of the orbital spin 
polarizations. 

To obtain further insight into the type of states associ- 
ated with each orbital Fig. [3] shows the projected density 
of states of the Mn d-orbitals for a full anisotropy antifer- 
romagnetic calculation. The large peaks in the majority 
up spin clustered between -8 and -5 eV indicate the po- 
sition of the more localized t^g orbitals. There is also a 
sharp peak for d x 2_ y 2 near -5 eV, but otherwise broader 
structure. This sharp peak is a result of the distortion 
of the in-plane O-Mn-0 angle to ~80° which results in 
partial charge, associated with d xy in the case of perfect 
octahedra, being attributed to d x 2_ y 2. This has been 
confirmed by calculating the projected density of states 
for a crystal structure relaxed with the constraint that 
all octahedral angles remain 90° in which case the sharp 
d x 2_ y 2 feature is not present. In contrast the structure 
of the d z 2 partial DOS is very broad and likely to be 
dominated by bonding with the apical oxygen. 



A. Effect of Anisotropic Interactions on the Band 

Gap 

The optimized crystal structure obtained applying full 
anisotropy was retained and the DOS with antiferromag- 
netic order computed within the effective U approach, 
resulting in a band gap of 0.4 eV, as compared with 0.8 
eV when full anisotropy is employed. Here, the use of 
the same crystal structure allows us to discriminate be- 
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FIG. 3. (Color online) Projected density of states of the Mn 
d-orbitals for a full anisotropy antiferromagnetic calculation 
from Wien2k. Postive is spin up and negative spin down. 
The inset illustrates the distortion of the approximate MnOe 
octahedron with an in-plane O-Mn-0 angle of ~80°. 



tween electronic and structural effects, so that we may 
confidently say that the application of the effective U 
approach leads to a reduction in band gap of 0.4 eV. 
This effect can be understood within a Kanamori-typc 
model^ of the interactions within the localized d-orbital 
manifold, allowing us to straightforwardly rationalize the 
key physics occurring in the ab initio calculations. The 
Hamiltonian for the d-states may be written: 



H = f + Ei 



(4) 



Here, Hq f incorporates the effect of the local environ- 
ment (crystal field) of the cation and Ej represents the 
DFT+U interaction term in equation [TJ Since the Mn 4+ 
ions occupy approximately octahedral sites, we will as- 
sume that the highest occupied d-state is t% g , where al- 
pha denotes the spin index, with an excitation across 
the band gap into an unoccupied state. Conse- 
quently, the d to d bandgap may be approximated as 
the difference in eigenvalues between these two states i.e. 
E g =e(e< g *)-e(t%). 

The DFT+U potential matrix elements with the fully 
localized limit double counting correction have previously 
been expressed by Ylvisaker et al. (see eq. 25 of refi^2): 



Av ma =-(£/- J) 



Av. 



aniso 
ma 



(5) 



The first term encompasses the isotropic terms and 
is identical in the effective U and full anisotropy ap- 
proaches. The second term Au™ so incorporates all of 
the effects of anisotropic interactions in DFT+U and will 
be zero in the effective U case, so that the difference 
in band gap in the effective U and full anisotropy ap- 
proaches must be largely attributable to the effect of this 
term alone. Below we detail the action of the latter term 
in our Kanamori Hamiltonian for /3-Mn02- 



The eigenvalues of the model Hamiltonian may be eval- 
uated using the Slater- Janak theorem as the derivative of 
the total energy with respect to the occupation number 
of the state e(i) = dE d ^ l ° t ■ For the effective U case: 



e{t% g ) = 2U ES 
e(e g l ) = 3U ES - 



*CF 



Gap Ef™ = C7 Eff + A 



CF 



Here Aqf lS the octahedral crystal field splitting, 
full anisotropy: 



(6) 
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(7) 



Considered together, the difference in bandgap be- 

Ef™ = (235/273) J, 
1.2 eV as in the 



tween the two methods is E^ nl 



amounting to 1.03 eV with a setting J 
ab initio calculations. However, we note that the band 
gap possesses a charge transfer character, evident from 
the dominant weight of oxygen p-states present in the up- 
per valence bands of Fig.[2jc). This implies that the ap- 
parent variation in band gap between the two approaches 
will be due to relative motion of the upper d-cigenvalue 
alone, amounting to a halving of the value above to 0.5 
eV prior to comparison with the ab initio calculations. 
This is broadly consistent with the variation in band gap 
of 0.4 eV that we see in the ab initio calculations. The 
small difference is likely due to the neglect of the effects 
of differences in the degree of Mn-0 orbital hybridization 
between the two approaches. This stems from the fact 
that the Acf parameter above should more realistically 
be considered as a ligand field rather than a crystal field 
splitting, in that it arises from a combination of rigid ion 
electrostatic and hybridization effects. 



B. Noncollinear Magnetism 

The full helical magnetic order with a pitch of (7/2)c 
first determined by Yoshimori^ is also considered here, 
where we assume that any effects due to the 4% deviation 
from this pitch refined by Regulski et al. are smalts. A 
supercell consisting of a seven unit cell expansion along 
the c-axis was therefore constructed, wherein the heli- 
cally ordered moments lie within the ab-plane. Defin- 
ing a propagation vector k = (0,0, k z = 2/7), the Mn 
sites positioned at (0,0,0) in each of the original unit cells 
(z„ = n) take moments: 

M (o,o,o) = M[sin(27rM„),cos(27rfc z z„),0] (8) 

where n indexes the units in the supercell. Meanwhile, 
the (|, h, h) positions (z„ = n + 1/2) in each unit cell 
take moments: 



M[— sin(27rfc z z„), — cos(2irk z z n ), 0] (9) 
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The total energy of this full helical order relative to FM 
ordering has been computed using the non-collinear mag- 
netism approach within the VASP code. Application of 
the effective U method leads to the stabilization of FM 
order by 65 meV per formula unit, while the inclusion of 
anisotropic interactions, with U and J as defined above, 
results in a stabilization of the full helical order by 10 
meV per formula unit. Inspection of the DOS for these 
non-collinear calculations gives band gaps that differ by 
at most 0.1 cV from those shown in Fig.[5]for the collinear 
calculations. This is due to the low energy scale for the 
non-collinearity which results in no qualitative difference 
to our conclusions regarding the band gap. 

Yoshimori^ has studied the magnetic structure of /3- 
M11O2 using a classical Heisenberg spin Hamiltonian in- 
cluding the first three Mn-Mn neighboring interactions 
Jinn, J2NN and J3NN as shown in Fig. [TJa). The he- 
lical magnetic order observed experimentally was found 
to be stable if: (1) Jinn, J2NN and J3NN are all posi- 
tive (favoring AFM order), (2) Jinn > Jinn and (3) 
J3NN is small or negative compared to J2NN accord- 
ing to the condition 2JinnJ%nn I J2.NN < To com- 
pare our DFT results to these conditions we have cho- 
sen to map the DFT energies onto a classical Heisenberg 
Hamiltonia n 33 ! 34 : 

H = E + 2j2 y £JZ j S:-S? (10) 

a (id) 

where a labels progressively more distant shells of Mn-Mn 
interactions. Couplings out to fourth-nearest-neighbor 
are explicitly considered here, but, as will be shown be- 
low, the fourth neighbor interactions are small enough 
that they may reasonably be neglected. The second sum- 
mation runs over distinct spin pairs occurring at each 
appropriate shell separation, a, where Sf is the three 
dimensional spin vector on site i and Eq is a constant 
term that can be thought of as the non-magnetic en- 
ergy of the system. For example, the energy of FM 
order within this Hamiltonian is obtained as Efai = 
E + [2J 1nn + 8J 2 nn + 4:J3NN + 8J4nn}S 2 . The Heisen- 
berg model was least-squares fitted to the total ener- 
gies of 13 different Ising-type spin configurations repre- 
sentable within a 2 x 2 x 2 supercell expansion of the ap- 
propriate relaxed crystal structure presented in Table [U 
To be clear, Ising-type configurations restrict the spin 
pairs occurring above to be fully co- or anti-aligned. The 
results from Wien2k and VASP for both effective U and 
full anisotropy are presented in Table IIIII The table in- 
cludes standard errors for all parameters obtained from 
the fit and the magnitude of the spin vector is S = 3/2, 
consistent with previous work and the d 3 formal configu- 
ration. Experimental results from the time-of-flight neu- 
tron study by Ohama et al. are included for comparison^ 
where a 15 % error was inferred from the fit to the ex- 
perimental data. 

Examining first the results from Wien2k, the effective 
U calculation yields AFM couplings for the first three 
nearest neighbor couplings. However, both Jinn and 



J2NN amount to less than half of the values obtained 
experimentally. Furthermore, while the effective U cal- 
culation satisfies Yoshimori's second condition for a he- 
lical state with Jinn / J2NN = 1-95 it grossly fails the 
third condition with a third nearest neighbor coupling 
that is too large resulting in 2J\nnJ3NN / J%nn = 4.26. 
Even within the range of error on this parameter, ±1.62, 
this places the spin system far from the stability condi- 
tion for the helical state. Meanwhile, using Wien2k AFM 
couplings are again obtained from the full anisotropy cal- 
culation for the first three neighbors, but are significantly 
larger in magnitude. The nearest two neighbor couplings 
deviate by less than 17% from those obtained experimen- 
tally and satisfy Yoshimori's second condition for a heli- 
cal state with Jinn / J2NN = 1.41. The third condition 
is borderline at 2J\nnJ3Nn/J'2nn = 1-05, but the large 
error renders it difficult to make a definitive statement 
as to the smallness of J3NN- Clearly, the inclusion of full 
anisotropy improves the description of the spin system 
by increasing the magnitude of J\nn and J2NN, as well 
as bringing the system to the border of the third stabil- 
ity condition for helical order. The corresponding mean 
field transition temperature may be calculated as TV = 

{S(S + l)/3k B ] {&Jl NN /JlNN + iJlNN ~ 8J 3 JViv) » 66 

K, in reasonable agreement with the experimental tran- 
sition of 92 K. Another important finding is that the 
fourth nearest neighbor coupling J^n n is small in all cal- 
culations, and always occurs with a large fractional error, 
suggesting that its contribution to the magnetic energy 
is ill defined. 

Turning to examine the VASP results, the effective U 
calculation with PAW potentials produces FM couplings 
for the first two nearest neighbors and therefore clearly 
fails Yoshimori's first and most basic condition determin- 
ing the stability of helical order. It is interesting that 
the effective U calculation using full-potential LAPW in 
Wien2k, by comparison, retains AFM couplings, albeit 
that they are small. This demonstrates that the approx- 
imations made to the potential in the PAW approach 
adversely affect the prediction of exchange constants in 
this system. A similar effect likely also underlies the dif- 
ference in magnitude of the energetic stabilities found 
for collinear order with effective U in VASP and Wicn2k. 
The full anisotropy calculation with VASP recovers AFM 
couplings for the three closest neighbors, although they 
are still significantly smaller than those found experimen- 
tally or those with full anisotropy in Wien2k, resulting 
in a low mean field transition temperature TV = 35 
K. Therefore, it seems clear that the inclusion of full 
anisotropy favors AFM coupling, but the apparent under- 
lying bias of the PAW potentials toward FM coupling re- 
sults in net smaller couplings in VASP as compared with 
Wien2k. The calculations in both Wicn2k and VASP 
are performed in the rotationally invariant scheme. The 
reasons for the poorer couplings from VASP warrants 
future work and may include not only the method for se- 
lecting the orbitals on which U acts, but also the details 
of the pseudization in the PAW potentials, the interplay 
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TABLE III. Exchange parameters in Kelvin obtained by fitting equation[lO]to total energies calculated for a range of magnetic 
states. Experimental results obtained by Ohama et al. are included for comparison^. Derived parameters Jinn / J2NN and 
2 Jijvjv JiNN / J2NN are n °t shown for the VASP effective U results on the basis that Jinn and J2NN do not satisfy Yoshimori's 
conditions for a stable helical stated. 





Wicn2k LAPW 


VASP PAW 


CRYSTAL09 


Experiment 




Effective U 


Full Anisotropy 


Effective U 


Full Anisotropy 


B3LYP 


PBEO 




JlNN (K) 


4.3±0.9 


9.0±1.7 


-19.1±6.2 


4.3±0.7 


13.2±2.1 


8.1±1.5 


8.9 


J2NN (K) 


2.2±0.3 


6.4±0.6 


-11.5±2.1 


3.6±0.3 


8.3±0.7 


5.8±0.5 


5.5 


J3NN (K) 


2.4±0.6 


2.4±1.2 


0.0±4.3 


1.8±0.6 


4.0±1.4 


3.2±1.0 


-1.3 


Jan n (K) 


-0.9±1.5 


-0.4±0.6 


-1.2±2.2 


-0.5±0.3 


-0.3±0.5 


-0.4±0.6 




Jinn/ J2NN 


1.95±0.48 


1.41±0.29 




1.19±0.21 


1.59±0.29 


1.40±0.28 


1.6^, 


2Jijvjv J3NN / J2NN 


4.26±1.62 


1.05±0.58 




1.19±0.46 


1.53±0.62 


1.54±0.59 


-0.76 



a The ratio 1.6 was assumed in fitting to the neutron time of flight experiments. 



between the Hubbard-like Hamiltonian and the pseudo- 
projectors and the use of a second variational procedure 
for the incorporation of the DFT+U potential. 

At this point it is timely to make some further general 
comments regarding the fitting of the exchange couplings. 
Firstly, theory in all cases predicts a positive (AFM) 
J3NN, while the experimental values yield a small neg- 
ative (FM) coupling J^nn = —1.3 K. However, it must 
be emphasised that the experimental results were ob- 
tained enforcing a ratio Jinn / J2NN = 1.6 in the fitting, 
which has an unknown effect on the J^nn value. Sec- 
ondly, four-parameter fits in which J^nn — is enforced 
showed slightly decreased values of Jinn and J3NN with 
little change in J2NN- For example, the full anisotropy 
calculation in Wien2k then gives Jinn = 8.4 ± 1.3 K, 
J2NN — 6.3 ± 0.5 K and J%nn = 1-9 ± 0.8 K. In sys- 
tems possessing rich magnetic phase diagrams, such as 
the present /3-Mn02, variations of this magnitude could 
conceivably alter the nature of the predicted low temper- 
ature ground state, in particular here whether the stabil- 
ity conditions set out by Yoshimori for helical order are 
satisfied. However, such small changes in coupling do 
not alter the overriding conclusion that fully anisotropic 
interactions improve the description of the exchange con- 
stants. 

It may also be suggested that the description of the 
magnetism and band gap could be improved by the 
use of a smaller value of U within the effective U ap- 
proach. Indeed the antiferromagnetic band gap only 
closes for values of Ues below 3 eV. Employing Ues = 3.5 
eV one may stabilize the collinear antiferromagnetic or- 
der by very small energies while still retaining a band 
gap. However, a full fitting of the exchange integrals in 
Wien2k with Ues = 3-5 eV results in Jinn — 8.5 ± 1.1, 
J 2NN = 4.0 ± 0.4 and J 3N n = 2.6 ± 0.8. With these val- 
ues J\nn / J2NN = 2.12 ± 0.37 and 2JinnJ3Nn/ J%nn = 
2.75±0.92 clearly demonstrating that Yoshimori's stabil- 
ity conditions are not satisfied. Therefore, reducing the 
value of Ues does not remedy the poor description of the 
exchange couplings. Furthermore, such a small value of 
U is not consistent with our self-consistent calculation of 
the parameter. 



Finally, for comparison, in Table Mil the results ob- 
tained from two types of hybrid functional calculations 
within the CRYSTAL09 LCAO code^ are also pre- 
sented. Here, the crystal structure obtained from hy- 
brid HSE06 3? in VASP is used giving a = 4.383 A and 
c = 2.875 A in good accord with experiment. The well- 
established 8-6411d41 and 6-31dl basis sets were em- 
ployed for Mn and O respectively^, other numerical pa- 
rameters being similar to the VASP and Wien2k calcula- 
tions. The results of both functionals satisfy Yoshimori's 
first and second condition for the helical order, although 
B3LYP (bearing 20% HF exchange)^ does overestimate 
the first three nearest neighbour couplings compared to 
experiment. PBEO (25% HF exchange)^ gives a similar 
level of agreement as was obtained with full anisotropy 
in Wicn2k. For both B3LYP and PBEO, the computed 
ratios 2J\nn Jznn/ J%nn = 1-53 and 1.54, respectively, 
appear to be too large to satisfy Yoshimori's third con- 
dition, but, again, the large errors hinder a definitive 
statement. 

The good performance of these hybrid functionals is 
consistent with the improvement obtained by use of fully 
anisotropic DFT+U over the effective U approach. This 
is not surprising, given that the proper anisotropy of d-d 
exchange interactions is included by construction in both 
hybrids, and that these are precisely the interactions de- 
scribed by full anisotropic DFT+U that do not occur in 
the effective U approach. Viewed more broadly, this sug- 
gests that, if the results obtained from hybrid functionals 
are to be compared with those from DFT+U, the latter 
calculations should apply fully anisotropic d-d interac- 
tions. 

We digress to discuss some qualitative features of the 
computed exchange couplings, particularly the antiferro- 
magnetic Jinn interaction. Antiferromagnetic coupling 
is predicted by all of our calculations (except effective 
U with PAW potentials) and is consistent with the ex- 
perimental findings. However, the Mn-O-Mn bond an- 
gle of the Jinn interaction, at ~101°, is close to 90°, 
and might therefore be predicted to yield weakly ferro- 
magnetic coupling based upon a straightforward appli- 
cation of the Goodcnough-Kanamori rule^ for 90° su- 
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perexchange. Beyond this simple model, the situation in 
the real crystal is made more complex by the further note 
of Goodcnough that "the cation-cation interactions tend 
to dominate the 90° superexchange in oxides if the t2 9 
orbitals are half-filled" . Thus, direct exchange should be 
considered in addition to the oxygen mediated superex- 
change, where the former may be ferromagnetic or anti- 
ferromagnetic depending upon the relative weight of the 
potential and kinetic contributions (varying sensitively 
with the orthogonality, separation and radius of the or- 
bitals involved). In light of the above, we suggest that 
an electron hopping driven antiferromagnetic direct ex- 
change dominates the Jinn coupling, as supported by 
additional calculations wherein the c-axis length of j3- 
M11O2 is reduced such that the Mn-O-Mn angle closes to 
90°. In that case Jinn becomes even more antiferromag- 
netic at 139 K using full anisotropy GGA+U, in keeping 
with an antiferromagnetic direct exchange that tends to 
be progressively strengthened by a reduction in cation 
separation. By way of contrast, the 180° Goodenough- 
Kanamori rules are based on lower order perturbation 
theory and might therefore be expected to be more ro- 
bust to changes in geometry, in keeping with our finding 
of antiferromagnetic J2NN couplings mediated by ^130° 
Mn-O-Mn bridges. 

Finally, Geertsma and Khomskii^ have also suggested 
that the position of a "pendant" ligand neigboring the 
bridging oxygen anions may induce antiferromagnetic 
90° coupling. They considered CuGeOs, where the Ge 
cations are pendant to the Cu-O-Cu bridges formed by 
Cu02 ribbons with nearly square planar anion coordi- 
nation of the Cu 2+ cations. Meanwhile, in the present 
/3-Mn02 phase, the Mn-O-Mn bridges arc formed by 
columns of MnOg polyhcdra edge sharing along the 
axis, while the pendant species comprise the Mn cations 
in neighboring chains. We have quantitatively tested the 
effect of varying the pendant Mn to bridging oxygen dis- 
tance while maintaining the bond lengths and angles of 
the Mn-O-Mn Jinn bridges at constant values. We find 
that a reduction in this distance of 0.1 A in fact de- 
creases Jinn in magnitude by 4 K using full anisotropy 
GGA+U, which is contrary to the trend expected if the 
pendant cations were promoting antiferromagnetic ex- 
change. This reinforces the suggestion above that the 
antiferromagnetic sign of Jinn is largely the result of di- 
rect Mn-Mn exchange. We stress again that the majority 
of our computed results are consistent with experimental 
findings for /3-Mn02- 

IV. CONCLUSIONS 

In summary, the primary outcome from this work is 
that the inclusion of fully anisotropic interactions in 
DFT+U is key to understanding the magnetic order and 



resulting band gap of /3-Mn02- This allows us to pin- 
point the origin of the low temperature magnetic tran- 
sition and resulting poor electrical conductivity that is 
observed experimentally. The isotropic component of 
d-d interactions is typically not treated adequately by 
GGA or LSDA, and there is no reason to expect that 
the treatment of anisotropic interactions should be any 
better. Therefore, the incorporation into DFT+U of fully 
anisotropic Coulomb and exchange interactions should be 
routinely considered for most d-clcctron systems, partic- 
ularly in light of the fact that the more accurate approach 
may be deployed with little additional computational ef- 
fort in the majority of electronic structure codes. More 
generally, these results are expected to be of importance 
at the interface of DFT+U methods with other treat- 
ments of correlated electrons beyond pure DFT, such as 
hybrid functionals, in which the anisotropy of such inter- 
actions is included by construction, and DFT+DMFT, 
where anisotropic interactions enter via a model Hamil- 
tonian. 

The mechanism leading to an increase in band gap due 
to the presence of anisotropic interactions has been ana- 
lytically demonstrated by considering the eigenvalues of 
a model Kanamori Hamiltonian. The Mn rf-orbital par- 
tial charges suggest that the improved treatment of full 
anisotropy is driven by a reduction in unphysical spin po- 
larization of predominantly bonding orbitals. Exchange 
couplings have been obtained by the fitting of a Heisen- 
berg spin Hamiltonian to the total energies of a range 
of magnetic states, and are found to be adversely af- 
fected by the pseudization in PAW potentials compared 
with full-potential LAPW and all-electron LCAO hybrid 
functional results. The exchange couplings obtained with 
DFT+U including full anisotropy and those from hybrid 
functionals produce the closest agreement with experi- 
mental results. The proximity of these calculated results 
to the threshold of the Yoshimori stability conditions for 
the helical order in /3-Mn02 suggest that this system is 
close to the border of competing magnetic instabilities. 
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